FFLO state in thin superconducting films 
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We present the analysis of the inhomogeneous Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) super- 
conducting state in thin superconducting films in the parallel magnetic field. For the tetragonal 
crystal symmetry (relevant to CeCoIns - the most probable candidate for the FFLO state forma- 
tion) we predict a very peculiar in-plane angular dependence of the FFLO critical field due to the 
orbital effect. In the uniform superconducting state the critical field should be isotropic. The mag- 
netic field pins also the direction of the FFLO modulation permitting thus to study the critical 
current anisotropy. Our calculations reveal a strong critical current anisotropy in the FFLO state 
in sharp contrast with the usual superconducting state. The predicted characteristic anisotropies 
of the critical field and critical current may provide an unambiguous probe of the FFLO phase 
formation. 

PACS numbers: 74.45.+C, 74.78.Fk 85.25. Cp, 74.81.-g 



I. INTRODUCTION 

Recently the strong experimental evidence has been obtained in favor of the existence of the inhomogeneous su- 
perconducting state in the heavy fermion compound CeCoIns (see ref. [l| and references cited therein). This state has 
been predicted a long time ago by Larkin and Ovchinniko\«2, and Fulde and Ferrell 3 - who demonstrated that in a pure 
ferromagnetic superconductor at low temperature the superconductivity may be non-uniform (for a review see refs. [l| 
and0 ). The conditions of such FFLO state formation are rather stringent. In CeCoIns the FFLO state exists owing 
to the paramagnetic pair-breaking effect which dominates the orbital one. Moreover the superconductivity of CeCoIn 5 
corresponds to the clean limit. Note also that the magnctic-ficld-induced superconductivity has been observed in the 
quasi two-dimensional organic conductor A-(BETS)2FeCl4, which may be another good candidate for the FFLO state 
realization^. 

Although the upper critical field in CeCoIns is dominated by the paramagnetic effect, the orbital effect also plays an 
important role in this compound leading to the vortex state formation. According to the theoretical works^iLSjSiifl in 
the FFLO state the vortex lattices may be very special and this circumstance complicates an unambiguous identifica- 
tion of the FFLO phase. The ideal system to study FFLO state without masking effect of the orbital field may be a two 
dimensional superconductor in a parallel magnetic field. It is a challenging experimental task to fabricate thin films of 
CeCoIns with thickness smaller than superconducting coherence length. There are strong experimental evidences that 
CeCoIns is a d-wave superconductor—. The potential scattering is harmful for both o?-wave superconductivity^ and 
FFLO phased. Therefore the films must have an atomically smooth surface and be in the clean limit. However due to 
the recent progress of the multilayer fabrication technology (in particular due to the epitaxial technique) such systems 
would be created in a near future. In the present work we study theoretically the properties of the FFLO state in 
superconducting film with the thickness d smaller than the superconducting coherence length £. To be more specific, 
we assume the CeCoIns multilayered structure with the tetragonal crystal symmetry and c axis perpendicular to the 
superconducting layers. Our approach is based on the modified Ginzburg-Landau (MGL) theoryii which adequately 
describes the FFLO state with a long wave length modulation - near the tricritical point (TCP). This TCP is the 
meeting point of three transition lines separating the normal, the uniform superconducting and the FFLO states. 
However qualitatively our results may be applied to the whole region of the FFLO state existence. The orbital effect 
for the magnetic field parallel to the film is weak, and comparing with the bulk superconductor its contribution to 
the upper critical field is reduced by the factor (d/£ c ) 2 <C 1, where £ c is the superconducting coherence length along 
c axis. Then we may expect a much larger FFLO region on the (H,T) phase diagram. In the present work, we 
demonstrate that small but finite orbital effect leads to a characteristic in-plane anisotropy of the upper critical field 
which may serve as a clear indication of the FFLO state formation. Also the magnetic field lifts the degeneracy of the 
direction of the FFLO modulation and permits to create a monodomain FFLO state. Moreover our analysis shows 
that we may expect a pronounced anisotropy of the superconducting critical current in FFLO state. 
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II. GENERALIZED GINZBURG-LANDAU FUNCTIONAL 

The long period FFLO modulation near the superconducting transition can be described by the MGL functional 
which in addition to the usual gradient terms contains the higher derivatives of the superconducting order parameter 
The necessity to add the higher derivatives is related to a special behavior of the coefficient on the usual gradient 
term which goes through zero and becomes negative in the region of FFLO phase. This circumstance makes the 
MGL theory qualitatively different from the standard Ginzburg-Landau approach and is responsible for the peculiar 
properties of the FFLO state. As an example, we note that in the Ginzburg-Landau functional the higher derivatives 
describe the weak non-local effects (see, for example, ref. [l3l ) which are related to the details of the Fermi surface. In 
the FFLO phase these non-local terms are of the primary importance and then making the properties of the FFLO 
phase ultimately dependent on the details of the electronic spectrum. 

The MGL functional provides an adequate description of the FFLO state near the trictitical point where the wave- 
vector of the FFLO modulation is small but obtained results could be extrapolated qualitatively on the whole region 
of the FFLO phase. In the case of pure paramagnetic effect this TCP corresponds to T* = 0.56T c o, H* — H(T*) = 
0.61Aq//!b = 1.05T c o//i_B, where T c0 is the critical temperature in the absence of the paramagnetic effect. The orbital 
effect decreases T* and in CeCoIn 5 T* seems to be (0.2 — 0.3)T c0 (ref. [l]) depending on the field orientation. In 
the bulk superconductors, the relative contribution of the paramagnetic and orbital effects may be characterized by 

/— H orb 

the Maki parameter a = V2 , where H°2 is the orbital critical field extrapolated to T = from the slope of 

p 

H C 2(T) near T c and H p is a paramagnetic limit at T = 0. In CcCoIns the Maki parameter is large a « 5 ensuring 
the condition for the FFLO phase formation. The particularity of CeCoIns is that the superconducting transition is 
slightly first order— below ~ 0.4T c o for the field H || ab and below ~ 0.3T c o for H _L ab, and it remains first order 
at FFLO transition^. This intriguing behavior may be related with rather pronounced magnetic fluctuations in this 
compound which is expected to be at the vicinity of quantum critical points. More generally the internal field h 
acting on the electron spins have the contribution coming from the Ce band, in addition to Zeeman's term JlH (where 
Jl is effective electron magnetic moment). If the exchange integral describing the interaction between this band and 
the superconducting electrons is /, the contribution from the polarized Cc atoms being SI, where S is their relative 
magnetization. Owing to the intcrband interaction, it would be a contribution to the Ce polarizability from the 

electron susceptibility x e i which changes at the superconducting transition \ e ~ X°, ~ c ^ ^A~J > where the constant 

c l ^ liS. In the result the internal exchange field h will acquire a correction ~ —IH \^\ . If the exchange integral 
is positive it will give a negative ~ |W| contribution in the Ginzburg-Landau functional. Note that for the negative 
exchange integral we may have an inverse situation and then the superconducting transition could be the second order 
one at all temperature. Therefore depending on the sign of the exchange integral, the Ce band could favor the first 
or second order superconducting transition. 

Taking CeCoIns in mind, we assume here that the superconducting transition is weakly first order and then the 
MGL functional reads: 

f g = a(i/,r)|*| 2 -a(|n^|V|n^| 2 ) +/3(|n^|V|n^| 2 + |n 2/ n ;i ,*| 2 + |n ;c n 2/ *| 2 ) (i) 
+£ (n 2 *(n 2 vi;)* + c .c.)-||*| 4 +^l*l 6 , 

where 11^ = ^ — i2eA x and Ii y = -g- — i2eA y , the superconducting film in xy plane and z axis perpendicular to the 
film. In the FFLO region the coefficients a, (3 > and the choice b, A > ensures the first order transition. Magnetic 
field is directed along the film and makes the angle 9 with x axis and then 

H y = Hsm9, H x = HcosO, 
A y = — zH cos 6, A x = zHsin.9. 

The coefficient a(H, T) vanishes at the line of the second order transition to the uniform superconducting state and 
at fixed H it may be written as a = ag(T — T CU (H)), where T CU {H) is the second order transition temperature into 
the uniform state. Due to the small thickness of the film d <C £ c , the superconducting order parameter is constant 
over its thickness and then in ([I]) we have omitted the derivatives on z. The higher derivatives terms with the 
coefficient e describe explicitly the difference between isotropic s-wave pairing model and the real situation realized 
in the tetragonal crystals and/or with d-wave pairing. In isotropic s-wave superconductor there is a degeneracy 
over the direction of the FFLO modulation. For d-wave superconductor this degeneracy is lifted and in 3D case 
the modulation vector is always directed along the order parameter nodes^i. In 2D d-wave superconductor at low 
temperature T < 0.06T C a first order re-orientational transition occurs to the state with a modulation vector along 
the order parameter lobesi^. Note that in general e ~ (3 and the effect of anisotropy for the FFLO state cannot be 
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expected to be small. This is an important difference with the standard Ginzburg-Landau theory where only the first 
rank tensors enter as a gradient terms and then the cubic crystal structure (or tetragonal in ab plane) is equivalent to 
the isotropic one. Then the form of the Fermi-surface and the type of the superconducting pairing are both equally 
important in determining the wave-vectors of FFLO modulation. For the FFLO transition the interplay between the 
Fermi surface structure and the type of the superconductivity has been studied in refs. and on the basis 
of the tight binding model and a very rich variety of the scenarios of the FFLO transition was revealed. In general 
it may be demonstrated^ that the effective mass approximation can be reduced to the isotropic model by scaling 
transformation. However there are namely the deviations from the elliptical Fermi-surface which are crucial to the 
adequate FFLO description. Note that the tensor coefficients on the second derivatives terms in MGL are given by 
the expression^ (3^ ~ (vfv 2 |-0(k)| 2 ), where Wj are the components of the Fermi velocity, ip(k) is the gap function 
and the averaging is performed over the Fermi surface. 



III. ORIENTATIONAL EFFECT OF THE IN-PLANE FIELD. IN-PLANE ANISOTROPY OF THE 

CRITICAL FIELD. 



Let us consider first the quadratic terms in ([T]) which depend on the orientation of the FFLO modulation. If 
the transition is of the second order (or weakly first order) the solution for the order parameter is of the form 
*( r ) = / cos(qr). 

Without orbital effect 

-aq 2 + f3q 4 + -q 4 sin 2 2tp 



8F{tp) 



(2) 



where tp is the angle between the FFLO modulation vector and x-axis. For e > the minimum energy (and maximum 
critical temperature) corresponds to tp = 0, ±7r/2, tt directions, i.e. along the x, y axis and the wave-vector go = \J~^- 

For e < the minimum energy corresponds to tp = ±37r/4, ±7r/4 directions, i.e. along the diagonals. 

In the case of a thin superconducting film in a parallel field it is easy to take into account the orbital effect - it is 
simply needed to average the functional ([T]) over the film thickness. The angular dependent part is 



SF(tp,6) 



+2q 2 (2eH) 2 ^ 
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(2eHY- + \ q * + {2eHy- 



so 



'- ( q 4 sin 2 2tp 



(2etf) 4 |jsm 2 20 



(3) 



2/3 + | - (j3 + |) cos(20 - 2tp) + ^ cos(20 + 2tp) 



Let us consider first the isotropic case e = 0. Naturally the properties of the superconducting system are not 
depending on the field orientation itself and the angular dependent part has the form — f3 \^ q \ q 2 d 2 (2eH) 2 cos(20 — 2<^) 
and the minimum of the energy is attended at = tp, i.e. when the FFLO modulation is directed along the magnetic 
field. So the magnetic field provides an orientational effect on the FFLO phase. In the absence of other sources of 
anisotropy the resulting field dependent contribution to the energy is isotropic in xy plane. The quadratic over H 

contribution vanishes for the wave- vector of the FFLO modulation q = qa = y^j- Then there is no linear diamagnetic 

response in the FFLO phase SF = [3 |* ? | 2 q 2 £■ (2eH) 2 - a(2eH) 2 |* g | 2 = . The resulting orbital field contribution 

~ \^ q \ 2 (3(2eH) 4 jgQ is quartic over H and then the diamagnetic moment is pretty small and proportional to H 3 . 

In the more realistic case the anisotropy (crystalline and/or Cooper pairing) plays a very important role and pins 
the orientation of the FFLO modulation. Let us suppose for example that e > and then in the absence of the orbital 
effect the FFLO modulation vectors are along the x, y axis. Note that in the case e < the rotation of the xy axis by 
45° provides us the same functional ((T|) with renormalized coefficients e and (3 but with e> 0. Therefore our analysis 
presented below may be directly adapted to this case. If (3 > | and the angle of the magnetic field \8\ < ir/4 then 
the direction of the wave vector q will be close to the x axis while for 37r/4 > \6\ > ir/A the modulation will be along 
y axis. For (3 < | the situation is inverse and the system chooses the modulation along the axis (x or y) making the 
largest angle with field direction. The deviation of the modulation direction from the principal axis x, y is small and 
for (3 > | (and \8\ < 7r/4) the equilibrium angle tp is 
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The diamagnctic moment of the FFLO state is strongly angular dependent. For (3 > | 

,2 OL 



M 







-Hd z 



/3(1 - |cos20|) + -(1 + |cos20|) 



(4) 



(5) 
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FIG. 1: Schematic presentation of the angular dependence of the in-plane critical field (or critical temperature) in FFLO state 
(solid line) and in the usual superconducting state (dotted line). 

Note that in the usual uniform superconducting phase there is no linear angular dependent contribution to the 
magnetic moment. 

The angular dependence of the critical temperature (critical field) for the FFLO state is 



a 



2 cr 



P + \-\{p-\) 



(6) 



This angular dependence is presented in Fig. 1. In the uniform phase the angular dependence appears only due to 



the anisotropy term and it is proportional to (2eH) : 

a(h,T)=a(2eHf^-P(2eH)^ 



e,„ „, 4 d 4 1 



{2eHf (1 - cos 4(9). 

2 80 2 



(7) 



Comparing the angular dependence of the critical field in FFLO phase (Eq.[S]with that in the uniform state (Eq.[7]), 
we see that the latter is much weaker (~ H 4 ) and has a different form - see Fig. 1. 

Therefore the experimental studies of the in plane anisotropy of the critical field above and below T* may provide 
a conclusive test of the FFLO state formation. As it has been already noted, the critical field in CeCoIns is mainly 
determined by the paramagnetic limit H ~ H p (the Maki parameter is large olm ~ 5). 

The orbital effect sufficiently far away from the tricritical point provides the relative contribution to critical field 

as cm (iTt) (io) ~ ao (a 1 )' A (|o) ' k ere a ° * s ^ e g ra dient term coefficient near T c q (i.e. far away from the FFLO 
transition) . The higher derivatives terms also contribute to the in-plane anisotropy of the critical field through the so 



of the order of 



where £o is the superconducting correlation length in plane. We may rewrite this as 



(&>) ~(^f(/o) 



with H. 



orb 



t- Near the tricritical point the usual orbital effect weakens 



called non-local corrections to the GL theory™ Their contribution is of the order of ( -rr 2 - i 

J \H orb J \£o J (a M ) \«o, 

Therefore the condition of the domination of the special FFLO behavior being — > , 1 ..> (4-) ■ Then even for 

1 b aa (a M ) \ Zo J 

d < the characteristic FFLO regime would be observed everywhere except a tiny vicinity of the tricritical point. 

In this section we supposed that the FFLO transition is a second order transition. In the case of the first order 
transition (like in CeCoIns ) the performed calculations give a field (or temperature) of the over cooling of the normal 
phase. The actual field of the first order transition will be somewhat higher. However in the case of a weak first order 
transition the corresponding field (critical temperature) is obtained by the simple shift of a(h,T) - see next section. 
Therefore in this case also we may expect a peculiar angular dependence (fig. I) of the critical field. 



(*) 
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IV. CURRENT IN THE FFLO STATE. ANISOTROPY OF THE CRITICAL CURRENT. 



As it has been discussed in the previous section the orbital effect of the parallel magnetic field is small but permits 
to orient the wave- vector of the FFLO modulation in combination with crystalline/pairing anisotropy and then the 
critical current of the film will be anisotropic. Near the transition into the FFLO state the minimum energy is achieved 
for the one dimensional cos-like modulation of the order parameter— i 10 ' 11 . Assuming j3 > § and the magnetic field 
oriented along the a;-axis we have 

*(r) = /cos(gox). (8) 

Apart from the choice of the direction of the FFLO modulation, the orbital effect of the parallel magnetic field 
leads only to the small rcnormalization of the coefficients of the ([1]). Therefore the FFLO in the thin film opens the 
possibility to study the critical current and its anisotropy. 

Naturally the ground state \&(r) = f cos(qox) has no current and to describe the current carrying states we choose 
the order parameter in the form 

*(r) = fcos(q x)exp(i<p(r)). (9) 

To calculate the in-plane current we need to introduce the parallel components of the vector-potential A|| = (A x , A y ) 
and the part of the functional describing the interaction with A\\ being 

25F A = a(H, T)f - af ( (Vy> - 2eA) 2 + g 2 ) + (10) 



-P ( 4 + 2e ^) 2 + 2 ?o (j% - 2*4,) 2 ) f + 2eql - 2eA y ) ' / 2 . 



In this expression we have retained only the leading gradient terms assuming 



<< qo and performed the av- 



eraging over the FFLO modulation. Note that for the current calculation we need to use its most general definition 
j = — and not the formula from the G 
namics of the MGL theory is in fact very 
the following expressions for the current 



j = — and not the formula from the Ginzburg-Landau theory. This is a consequence of the fact that the electrody- 
namics of the MGL theory is in fact very different from the standard GL ono22. Using the relation q% = ^ we obtain 



For ip = kx the state with the uniform current along the x axis is realized. It is interesting to note that the 
perpendicular current along the y axis (with (p = ky ) is proportional to the anisotropy parameter e and vanishes in 
the idealized isotropic model. 

Now we calculate the critical current in the FFLO state following the approach similar to the standard GL theory 
(see, for example, ref. H3) but taking into account the first order character of the FFLO transition. In the absence of 
the current the order parameter has the form ([5]) and putting it into ([T]) we obtain the averaged free energy density 

F G = (a-a 2 )(H,T)f- b -f + ±f, 

2 

where a = a 2 = corresponds to the second order transition into FFLO state. The first order transition occurs 
at higher temperature/magnetic field and its "temperature" a\ = a 2 + f§j and amplitude of the order parameter 
fo = fx ma y be easily found from the conditions 

F G (oi,/o) = (ai-a2)/ 2 -^/o 4 + ^/o =0, (12) 
9F G {ai,f ) 2 4 

TT72 = ("i - a V - bf + A/ = 0. (13) 

Note that the "critical temperature" a\ of the fist order transition is simply obtained from the "critical temperature" 
of the second order transition a 2 by the shift on ||j . 

It is convenient to introduce the normalized "temperature" and order parameter: t — {a — a 2 ) / (ai — a 2 ) , / = ///o- 
Without the current the "temperature" dependence of the order parameter is given by ff = (2 + V4 — 3t) /3. In a 
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FIG. 2: Superconducting current versus the normalized amplitude of the order parameter z = f/ft. Maximum current 
corresponds to the critical current of the FFLO state. For both orientations the maximum of the current is achieved at the 
same value z m 0.81. 



presence of the current along the x axis, taking into account the equations (p~0|) . (fTTj) . and (p~3|) we have the following 
relationship between the current and the amplitude of the order parameter. 

j 2 x = 3 -^f?ftz*(l-z*)(7-z% (14) 

where z = f/ft- The condition of the applicability of our approach \/tp << qo reads b 2 << aA/y 7 /?, i.e. the FFLO 

transition must be weakly first order. The plot j x (z) (sec Fig. 2) have a maximum at z 2 = (3^2- VTl) /v^^ 0.65 
which gives us the critical current along the x axis 

jL .„ ra H5^ /f/o , (U1) 

Analogously we find that the critical current in the direction perpendicular to the FFLO modulation is 

a _ 1 fe\ 2 2.85ae 2 b 2 T8H 

Jycrit — 4 \ /5 A JtJO- I 10 / 

The coefficient 1/4 is coming from the averaging on the order parameter modulation along the x axis. Finally the 
anisotropy of the critical current in the FFLO phase is very pronounced 

Jxcrit 

2(3 

Jycrit & 

and in the isotropic model the critical current along the y is even vanishing. In real system the ratio ^ is expected to 
be of the order of one and the measurements of the critical current can permit to determine directly this parameter. 
As the critical current in the uniform state is isotropic far away from the trictitical point, then the experimental 
observation of the anisotropy of the critical current may serve as a clear indication of the FFLO phase formation. 



V. CONCLUSIONS 



To summarize, we have investigated the properties of the FFLO phase in a thin film at parallel magnetic field. 
The orbital effect (even though it is small) leads to the orientation of the FFLO modulation through the whole film 
providing the monodomain FFLO state. Moreover in the FFLO state a peculiar angular dependence of the in-plane 
upper critical field must be observed. This conclusion is quite general and holds for both first and second order FFLO 
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transitions. We predict also an important anisotropy of the in-plane critical current in the FFLO phase depending 
on the current direction with respect to the FFLO modulation. Such characteristic anisotropies of the critical field 
and critical current may be considered as a smoking gun of the FFLO phase formation. Note that ir*2^ it has been 
demonstrated that the superconducting fluctuational regime changes drastically near the FFLO TCP. However in the 
case of the first order FFLO transition the fluctuational regime could be inaccessible on experiment. Our analysis 
was based on the very general MGL functional approach which is valid for both s-wave and d-wave superconductors. 
This approach is fully justified near the tricritical point and in the case of the weakly first order phase transition. 
Nevertheless qualitatively the obtained results could be extrapolated to the whole region of the FFLO phase existence 
and may be relevant for the CeCoIns thin film experiments. In conclusion we stress that the predicted anisotropy of 
the of the in plane critical field and critical current must be also observed in s-wave tetragonal FFLO superconductor. 
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